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ABSTRACT 

We measure the evolution of the rest-frame i^-band Fundamental Plane from z = 1 to the present 
by using IRAC imaging of a sample of early-type galaxies in the Chandra Deep Field-South at z ^ 1 
with accurately measured dynamical masses. We find that M/Lk evolves as A In [M/Lk) = (— 1.18± 
0.10)z, which is slower than in the B-band (Aln(M/Li3) = (-1.46 ± 0.09)z). In the B-band the 
evolution has been demonstrated to be strongly mass dependent. In the iC-band we find a weaker 
trend: galaxies more massive than M = 2 x lO^^Af© evolve as Aln (M/L/^) = (—1.01 ± 0.16)z; 
less massive galaxies evolve as A In {M/Lk) — (—1.27 ± 0.11)z. As expected from stellar population 
models the evolution in M/Lk is slower than the evolution in M/Lb- However, when we make a 
quantitative comparison, we find that the single burst Bruzual-Charlot models do not fit the results 
well, unless large dust opacities are allowed at z = 1. Models with a fiat IMF fit better, Maraston 
models with a different treatment of AGB stars fit best. These results show that the interpretation 
of rest-frame near-IR photometry is severely hampered by model uncertainties and therefore that the 
determination of galaxy masses from rest-frame near-IR photometry may be harder than was thought 
before. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation 



1. INTRODUCTION 

The formation and evolution of early-type galax- 
ies (hereafter, E/SOs) has been a major subject of 
research during the past decades. Both theoreti- 
cal and empirical studies indicate that the stellar 
mass density of the ETG population has increased 
by at least a factor of two between z = 1 and 
the present day fKauff mann fc Charlotill998t iBell et all 
|2Q04; Fabcr et al. 200^ At the same time, studies at 
z ~ 1 show that massi ve E/SOs were already several 
Gyr old at that epoch llvan Dokkum fc Stanfordll2?M 
Holden et aO 120051 llVeu et alJl2005t Ivan der Wei et all 
2005^ . It is an intriguing puzzle to explain these seem- 
ingly contradictory results in one picture describing 
early-type galaxy formation. 

These results hinge on determining galaxy masses, 
which remains very hard to do at high redshift. Typ- 
ically, mass estimates for high-z galaxies are obtained 
by comparing their photometric properties with stellar 
population models in order to infer a mass-to-light 
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ratio (M/L) (see, e g.. IBell et a1.l HOfll IBuudv et a,1J 
l2005t IShaolev et all mm) . Such SED-fitting meth- 
ods are intrinsically uncertain, and dynamical mass 
calibrations are needed to overcome the large un- 
certainty in the low mass end of the initial mass 
function (IMF), dark matter content, and, most im- 
portantly, to verify the validity of stellar population 
models. A more quantitati ve approach is to use the 
Fundamental Plan e (FP, iDiorgovski fc Davis! llQSTt 
iDresslcr ct al.' "1987') to directly measure the evolution 
of M/L at rcst-framc optical wave lengt hs out to z = 1.3 
jvan Dokkum fc Franx .19961 Iva n Dokku m et al.l llQQSt 
Kelson et aL.200q: .Treu et al.ll200L .van Dokkum et all 
200n iTreu et a,].l 120021 lva,n Dokkum fc Sta.nfordI [200? 
'van Dokkum fc Ellis' '2003; 'Gebhardt et al.1 120^ 
van dcr Wcl ct al. 2004; Wuvts ct al. 2004; Hol den et all 
|2005; Treu ct al. 2005; van der Wei et al.ll200^ . 

Although the measurement of the evolution of M/L it- 
self is model independent, the inferred age and formation 
redshift are sensitive to the choice of the parameters (e.g., 
IMF and mctallicity) of the stellar population model, and 
even the choice of the model. By measuring the evolu- 
tion of M/L at other wavelengths than the rest-frame B- 
band, the formation redshift can be established indepen- 
dently, and the range of possible model parameters con- 
strained. Furthermore, the validity of different models 
can be verified. More specifically, different models pre- 
dict a different evolution of the M/L in the near infra-red 
(NIR) relative to the evolution of the optical M/L. An 
advantage of measuring the evolution of the NIR M/L is 
that the NIR luminosity is less affected by the presence of 
a low mass young stellar population. Since measuring the 
evolution of M/L provides a luminosity- weighted age es- 
timate, the evolution of the NIR M/L provides a less bi- 
ased age estimate. Wit h the arrival of the near Infra-Red 
Array Camera (IRAC, iFazio et aDl2004() on the Spitzer 
Space Telescope, studying the rest-frame NIR proper- 
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ties of the high-z E/SO population has become feasible. 
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Fig. 1. — SEDs of three z ~ 1 E/SOs. The red asterisks are 
observed fluxes in respectively the ACS («, i, z), ISAAC (J, 
K), and IRAC (3.6, 4.5, 5.8, 8.0) bands. The lines are model 
spectra from Bruzual- Chariot (a 2 Gyr-old, dust-free SSP 
with solar metallicity and a Salpeter IMF) The photometric 
SEDs, including the IRAC data points, have plausible shapes 
of stellar populations of several Gyrs old. 

In this Le tter we combine IRAC imaging with our high-z 
FP study jvan der Wei et alJl200ll2005l) . such that the 
evolution of M/L in the rest-frame fsT-band is measured 
for the first time. A calibration of galaxy masses derived 
from NIR luminosities is thus provided. 



2. IRAC PHOTOMETRY OF E/SOS AT Z ~ 1 

Ivan der Wei et all l)2005D provide a spectroscopic sam- 
ple of 29 distant galaxies in the CDFS with accurate 
velocity dispersions. Spitzer GTO data are available, 
providing IRAC imaging in four channels (3.6/i, 4.5/i, 
5.8/i, and 8.0/i) of a field containing the Chandra Deep 
Field-South (CDFS). Using the IRAC data we derive the 
rest-frame NIR photometry of those 20 galaxies in the 
spectroscopic sample with early-type morphologies and 
sufficiently high S/N spectra (S/N > 12 per A). The 
average redshift of this sample is z = 0.94 and the aver- 
age mass is 10^^ Mq. 12 galaxies in this sample are re- 
garded as the 'primary sample': these objects are all at 
0.9 5 < z < 1.15 and s atisfy all selection criteria applied 
bv lvan der Wei et al.l (j2005). The other objects, mostly 
galaxies at z 0.7, are referred to as the 'secondary 
sample'. The IRAC data are sufficiently deep (500s) to 
obtain high S/N photometry of all galaxies in the sam- 
ple, certainly in the two shortest wavelength channels. 

The IRAC fluxes were measured by matching the 
point spread functions of all available optical/NIR/IRAC 
images to the S.Ofi IRAC image (which has the low- 
est spatial resolution, FWHM = 2'.'3). 5'.'0 diam- 
eter aperture fluxes we re measured using SExtractor 
ijBertin fc Arnout^ Il99()j) to derive colors. Color gra- 
dients are ignored, but because the sizes of the z ~ 1 
galaxies are much smaller (re// ~ 0'.'4) than the size 
of the IRAC PSF this will only have a very small ef- 
fect. The median formal errors in the IRAC fluxes 
of the galaxies in the sample are 0.01, 0.01, 0.10 and 
0.13 mag for the four channels, respectively. Addi- 
tionally, the systematic uncertainty in the IRAC zero- 
points is 0.03 mag. We transformed the observed col- 
ors to rest-frame B — K colors using a procedure simi - 
lar to what is described bv lvan Dokkum fc Franxl l)1996j) . 
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Fig. 2. — M/Lk evolution with redshift of E/SOs in the 
CDFS. The solid dots are objects in the primary sample, the 
crosses are the secondary sample. The red symbols are galax- 
ies with high masses (M > 2 x 10^^ M©), the blue symbols are 
less massive galaxies. The black data point labeled with 'Lo- 
cal' is derived from local cluster galaxies, but a small differ- 
ence between field and cluster galaxies is taken into account. 
The black line indicates the evolution of the entire sample, 
the red line indicates the evolution of massive galaxies in the 
primary sample. The evolution of M/ Lk of massive galax- 
ies is ^ 30% slower than the evolution of M/Lb of the same 
sub-sample. 

The typical error in the rest-frame B — K color is 0.06 
mag, excluding the systematic error of 0.03 mag. Total 
luminosities were derived from the structural parameters 
as obtained from the ACS z-band images and the colors 
are described above. 

From the higher resolution K-h&nd image we checked 
whether the IRAC photometry of our sample suffers from 
confusion. Three of the objects in our sample have 
close neighbors with such magnitudes that our photom- 
etry might be erroneous. However, we do not see a 
difference between these contaminated objects and the 
rest of the sample in either color, M/L, or any other 
parameter. In Figure 1 we show the SEDs of three 
high-redshift E/SO s . As an illustration we overplot a 
iBruzual & Charloil l)2003() model spectrum for a 2 Gyr- 
old stellar population, demonstrating that the observed 
SEDs have reasonable shapes. 

We also measured the colors of a sample of 41 E/SOs in 
the CDFS at 0.6 < z < 0.8. These were morphologically 
selected by eye from the ACS imagiiig. Th e redshifts are 
taken from COMBO- 17 l|Wolf et all 1200^ . We use the 
average B — K color of this sample to complement the 
B — K evolution of the galaxies in the FP sample. 

3. EVOLUTION OF M/Lk 

If galaxies evolve passively, the offset of high-z galaxies 
from the local FP is a measure of the difference between 
the M/L of the dista nt galaxy and the M/L and their 
local counterparts (see lvan Dokkum fc Franxll996j) . The 
rate of evolution is a measure of the relative age difference 
between the local and distant galaxies. 

To obtain the M/ L evolution in the iiT-band, we com- 
pare the NIR properties of the distant gal axies with the 
local FP in the if-band l|Pahre et al.lll99a() . i.e., we need 
effective radii [veff) and surface brightnesses {fieff) of 
the distant galaxies in the rest-frame iiT-band. How- 
ever, the spatial resolution of IRAC is too low to mea- 
sure galaxy sizes at z ~ 1. Instead we assume that 
feff,K = i^eff,B and determine the surface brightness 
at this radius. Although this causes some error in the 
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FP, previous work has shown that the combination of 
ref / and Ue Jj which enters the FP is ver y stable against 
such errors l)van Dokkum &: FranxllT996|) . The NIR sur- 
face brig htness within the o ptical effective radius is pro- 
vided bv lPahre et a l.' fl99^ for a local sample of cluster 
galaxies. For the distant sample we compute fieff,K from 
tJ-eff,B and (B — K)eff, where {B — K)eff is the color 
inside the effective radius. It is calculated from the total 
B — K color by c orrecting for the neg ative color gradi- 
ent measured by iPeletier et al.l 1)1990(1 . The difference 
between B — K and {B — K)eff is -0.07 mag. 

The local if-band FP is based on cluster galaxies, but 
our sample consists of field galaxies. iFaber et alJ l(1989l ^ 
have shown that field galaxies have ~ 5% lower M/Lb 
than cluster galaxies. Assuming that this difference is 
caused by either a difference in age or metallicity, this 
translates into a 2% difference in M /Lk, which follows 
from s tellar population models from lBruzual fc Chariot) 
li200l . We take this differ ence into account when com- 
paring the FP at high and low redshifts. 

We show the evolution of M/Lk with redshift in Fig- 
ure 2. It is apparent that M/Lk evolves significantly 
from z ^ 1 to the present. The evolution of the sam- 
ple is Aln(Af/iK) = (-1.26 ± 0.15)z, obtained from 
a least squares linear fit. The primary sample alone 
evolves at a similar rate: /S.\iv [M / L k) = (— 1.18±0.10)z. 
Comparing these numbers to the evolution of M/Lb 
{A\n{M/LB) = (-1.52 ± 0.16)z and (-1.46 ± 0.09)z, 
respectively), we see that the evolution in M/Lk is 
somewhat slower than the evolution in M/Lb- The 
scatter in Ahi {M/Lk)/z is 0.32, which is not much 
smaller than the scatter in the i?-band (0.37). This 
is much larger than can be accounted for by mea- 
surement errors, and is most likely caused by an age 
spread of the stellar populations of the galaxies, as has 
been demonstrated before (van dcr Wcl ct al. 2005). In 
the B-band the measured evolution is strongly mass- 
dependent llvan der Wei et all 120041: iTreu et all 120051 
Ivan der Wei et al.ll2005l^ . mostly due to selection effects, 
but partially due to intrinsic differences between high- 
and low-mass galaxies (van dcr Wcl ct al. 2005). In the 
K-hand the difference is much less pronounced, and only 
marginally significant. The galaxies in the primary sam- 
ple with masses higher than M — 2 x 10^^ Mq evolve 
as Aln(M/Li<-) = (—1.01 ± 0.16)2:; galaxies less massive 
than that as Aln(M/L/f) = (-1.27 ± 0.11)z. In the 
i?-band these numbers are, respectively, (—1.20 ± 0.14)2: 
and (—1.60 ± 0.09)2. These numbers are consistent with 
the model prediction that M/Lk is less sensitive than 
M/Lb to the age of a stellar population and recent star 
formation. 

Concluding, we find that the evolution of M/Lk from 
2 = 1 to the present is ~ 30% slower than the evolution 
of M/Lb- We note that this direct measurement de- 
viates from previous determinations of the evolution of 
M/Lk that were based on extrapolating observed SEDs 
(including the X-band as the longest wavelength data) 
to the rest-frame i^'-band at 2 1 using stellar popu- 
lation models. Such studies found that M /Lk evolves 
at lea st twice as slow as M/Lb (see, e.g.. iDrorv et al.l 
1200 4*). This demonstrates that extrapolating observed 
X-band photometry of high-2 galaxies to rest-frame K- 
band M/L is hazardous and strongly model dependent. 



4. DISCUSSION 

Our results allow us to compare the evolution of M/Lb 
and M/Lk directly with predictions from stellar popu- 
lation models. We compar e our results with the prc- 
dictions of three mo dels: iBruzual fc Chariot (^003), 
iVazdekis etHI (|1996|) . and iMarastonI l|2004D Chereafter. 
BC03, V96, and M04, respectively). A critical aspect 
of the models is the method that is used to implement 
late stellar evolutionary phases. BC03 and V96 compute 
isochrones up to the early AGB phase, and include an 
empirical prescription for the thermally pulsating- (TP- 
) AGB phase. On the contrary, M04 adopt the 'fuel 
consumption' approach, which allows implementation of 
short-duration but very luminous evolutionary stages, 
such as the TP-AGB phases, in an analytical and nu- 
merically stable way. 

In Figure 3 we compare the evolution from 2 = 1 to the 
present of M/Lb and M/Lk, and the related change in 
the rest-frame B — K color. The BC03 model with solar 
metallicity and a Salpeter IMF does not fit the data: the 
predicted evolution of M/Lk against M/Lb is signifi- 
cantly slower than observed. Before investigating differ- 
ent models and model parameters, we consider several 
possible explanations for an apparently fast evolution 
of M/Lk with respect to M/Lb that are unrelated to 
stellar populations. First, a difference in dust content 
between local E/SOs and 2 ~ 1 E/SOs could lead to red- 
der colors at 2 = 1 than expected from dust-free models. 
However, the amount of reddening required to match the 
BC03 model with solar metallicity and a Salpeter IMF is 
considerably large for half of the galaxies in our sample 
{Ay > 0.5), and several even require Ay > 1. Although 
the optical spectra show no indication of such high opaci- 
ties, observations at longer wavelengths with Spitzer can 
constrain the presence of dust-enshrouded populations 
which might contribute in the K-band. One might ex- 
pect to observe irregular morphologies or irregular inten- 
sity profiles if the absorption were so high. We note that 
the galaxies do not show such effects at the ACS resolu- 
tion. Second, the red colors might be produced by AGN 
in the 2 ~ 1 sample that are obscured in the optical, but 
not entirely in the NIR. 

Next, we test whether changing parameters within the 
single-burst BC03 model can improve the consistency be- 
tween the model and the data. Increasing the metallicity 
does not change the predictions by much. Changing the 
star formation history (SFH) does not lead to better fits 
to the evolution of M/Lb and M/Lk, unless models with 
constant star formation rates or exponentially declining 
star formation rates with e-folding time scales of several 
Gyr are adopted. These models, however, have very blue 
colors that do not match the B — K colors of the galax- 
ies in our sample, nor the colors of the local population. 
Flattening the slope of the IMF (from 1.35 to 0.35) does 
shift the model curve toward our data points, leading to 
an acceptable fit. 

Finally, we compare our results to other models. As 
can be seen in Figure 3, the V96 and M04 models with 
a Salpeter IMF and solar metallicity provide better fits 
than BC03. Note that the differences between the models 
with identical parameters are very large; e.g., BC03 pre- 
dict a correlation between color and A In (M/L), whereas 
M04 does not. We note that the spread of the data 
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Fig. 3. — Comparison of tlie M/L evolution in the rest-frame B- and i('-bands and the B — K colors of our z ~ 1 E/SO 
sample. (Left): the evolution of M/Lb versus the evolution of M/Lk, normalized at z = 1. The large filled symbol with 
error bars is the average of the four galaxies in the primary sample with masses larger than M = 2 x 10^^ Mq. The large open 
symbol with error bars represents the eight less massive galaxies in the primary sample. The smaller symbols are the individual 
galaxies in the primary sample. The diagonal error bar at the right shows the typical errors in AlnM/L_B and AlnAf/Lx for 
the individual data points. The shorter side of the error bar represents the uncertainty in B ~ K. The vector indicates the 
effect of dust. (Right): the evolution oi B — K versus the evolution of M/ Lk, normalized at z = 1. The thick horizontal error 
bar at the bottom indicates the average B ~ K evolution extrapolated to z = 1 of a sample of 41 E/SOs at 0.6 < z < 0.8. The 
model tracks indicate the expected evolution between z = 1 and the present, varying the formation redshift along the tracks. 
The Bruzual- Chariot models with a Salpeter IMF fit badly to the data. A Bruzual-Charlot model with a fiat IMF provides 
a better fit. The Maraston models with a Salpeter IMF provides the best fit. These models have a different treatment of the 
AGE stars. 



points representing the individual galaxies is large, but 
that there is no clear correlation between Aln (M/Lk) 
and A{B - K). Even assuming a flat IMF, the BC03 
model cannot reproduce the data points at the bottom 
right part of Figure 3b. The M04 model docs fit to those 
data points. On the other hand, data points at the upper 
left of Figure 3b are not fit well by the M04 model. How- 
ever, all but one of the these data points can be fit with 
a model with an exponentially declining star formation 
rate with r = — 1 Gyr. 

We can use the simple single burst models to esti- 
mate the average, luminosity-weighted formation redshift 
of the stellar populations of the galaxies. The forma- 
tion redshift as estimated from the evolution of M/Lb 
is 1.6 < Zform < 1-9 and 1.5 < Zform < 1-9, us- 
ing the models from BC03 and M04, respectively, with 
a Salpeter IMF and solar metallicity. According to 
the same models, the evolution of M/Lk suggests that 
1.2 < Zform < 1-5 and 1.6 < Zform < 2.0, respectively. 
The BC03 model with a flat IMF provides a better fit 
to our results than the BC03 model with a Salpeter IMF 
and produces 1.7 < Zform < 2.7, both from the evolu- 
tion in the M/Lb and the evolution in M/Lk- As can 
be seen, the constraints on Zform are not improved by 
measuring the evolution of M/Lk- This is entirely due 



to the rather large model uncertainty in the rest-frame 
K-hsiad evolution. Obviously, it is of the greatest rele- 
vance to improve our knowledge of both the AGB phase 
and the dust content of high redshift elliptical galaxies, 
as dust may mimic some of the results we see. 

Finally, we note that E/SOs are thought to be rela- 
tively simple systems, with simple SFHs. It is likely that 
the model uncertainties are even higher for galaxies with 
a more complex SFH. Irrespective of the causes of the 
discrepancies between the models and the observations, 
the results presented here suggest that the rest-frame 
K-hmid photometry of such galaxies may be difficult to 
interpret. Further studies of the role of dust and AGN 
would be valuable. If they do play a role in the pho- 
tometric bands analyzed here, proper modeling of the 
colors is more complex than thought before. In the next 
paper, we will explore full SED fitting, the uncertainties 
in determining stellar masses, and distinguish between 
the applicability of different stellar population models. 
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